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Abstract 

The peptide hormone cholecystokinin (CCK) exerts a wide range of digestive 
and CNS-related physiological signaling via CCK receptors in brain and gut. 
There is very limited information available on these receptors in Atlantic 
salmon. The aim of this study was to characterize CCK receptors in gut and 
brain of salmon. We have identified and cloned one CCK-1 receptor and 
duplicates of CCK-2 receptor in salmon. The phylogenetic analysis indicates the 
existence of one common ancestor gene for all CCK receptors. CCK-1R mRNA 
is highly expressed in pancreas followed by midgut, hindgut, gallbladder, and 
stomach indicating an involvement in pancreatic regulation and gallbladder 
contractions. CCK-2Rl/gastrin mRNA is expressed at high levels in midgut and 
at relatively low levels in stomach, gallbladder, and pancreas. We postulate 
CCK-2Rl/gastrin receptor to have gastrin-related functions because of its distri- 
bution and abundance in gastro-intestinal (GI) tissues. CCK-2R2 is relatively 
abundant in brain but has low expression levels in gut tissues supporting the 
hypothesis for involvement in the gut-brain signaling. Major functional motifs 
and ligand interaction sites in salmon are conserved with that of mammals. This 
information will be instrumental for comparative studies and further targeting 
receptor activation and selectivity of biological responses of CCK in salmon. 
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Introduction 

Cholecystokinin (CCK) is a key regulatory peptide 
hormone that serves a wide variety of physiological 
actions. It is a prototype of a class of agents known as 
brain-gut peptides, functioning both as neuropeptide and 
gut hormone. In peripheral tissues, CCK is known to be 
involved in multiple biological processes such as pancre- 
atic exocrine secretion, insulin release, and gut motility 
(Wank 1995; Tullio et al. 2000; Dufresne et al. 2006; 
Rehfeld et al. 2007; Stab'anssens et al. 2011). CCK also 
affects cell growth (pancreas and various tumors). Within 
the brain, CCK actions elicit satiation, anxiety-related 
behavioral responses and stimulate growth proliferation 



signaling protein extracellular signal-regulated kinase 1/2: 
ERK1/2 (Yassin 1999; Piiper et al. 2003). 

CCK exerts its biological actions by binding to specific 
G-protein-coupled receptors (GPCRs) known as CCK-1R 
and CCK-2R (formerly also termed as CCK-AR and CCK- 
BR, respectively). These receptors differ in their distribu- 
tion and their affinities for natural agonists. CCK-1R is 
more predominant in the periphery especially in pancreas, 
gallbladder, and gastric mucosa whereas CCK-2R is more 
abundant in different parts of central nervous system 
(CNS) (Tullio et al. 2000). CCK-1R is highly selective for 
sulfated CCK analogues. Gastrin receptor in stomach is clo- 
sely related to the CCK-2R subtype and they are often 
referred as combined CCKB/gastrin receptor as they show 
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similar high affinity for both sulfated and nonsulfated pep- 
tide analogues (Lee et al. 1993; Dufresne et al. 2006). 

For the CCK receptors in alimentary tract, comprehen- 
sive knowledge is available on their involvement in bio- 
logical processes in mammals, while there is very limited 
information on CCK receptors in fishes (Furutani et al. 
2013). The aims of this study were therefore to identify 
CCK receptors in Atlantic salmon and to elucidate their 
expression patterns in different gastro-intestinal (GI) 
compartments as well as in brain. In addition, we have 
assessed the putative phylogenetic relationships existing 
among CCK-R genes and predicted the 3D structure, 
interacting motifs and the agonist-binding sites of these 
receptors. Altogether these studies provide new insights 
into the conservation of CCK-R genes structure and most 
likely functions during evolution, as well as a better 
understanding for further studies in receptor activation 
and downstream biological responses of CCK. 

Material and Methods 

Specimen and tissues 

Atlantic salmon (average body weight 165 ± 3.97 g) were 
reared at the Bergen High Technology Center (Norway), in 
indoor tanks (500 L rearing volume) supplied with continu- 
ous water flow (2.5 L min -1 , 8°C) and fed ad libitum on a 
commercial diet (EWOS, Norway) for 3 weeks, as described 
by Ronnestad et al. (2010). All handling of animals followed 
the Norwegian Animal Research Authority (NARA) guide- 
lines for experimental purposes. At the end of the experiment 
(9 am, after feeding), fish were randomly sampled and sacri- 
ficed by an overdose of MS-222 (3-aminobenzoic acid ethyl 
ester; Sigma, St. Louis, MO). Tissue samples were collected 
in liquid nitrogen and kept at -80°C until further analysis. 

Gene expression studies 
RNA extraction 

Total RNA was extracted using TRI reagent (Sigma, St. 
Louis, MO). Isolated RNA was further purified by DNAse 
treatment (Turbo DNA free kit, Ambion, Life Technolo- 
gies Corporation, CA). First strand cDNA was synthesized 
from 4 /(g of total RNA using oligo(dT)20 primer and 
Superscript III (Invitrogen, Carlsbad, CA) according to 
manufacturer's protocol. 

Identification and cloning of Salmo salar CCK-1R, 
CCK-2R1, and CCK-2R2 genes 

Zebrafish CCK-1R and CCK-2R open reading frames 
(ORFs) (acc.no. NM_001085808.1 and XM-002663315.2) 



were used as a query against the codgenome project 
(http://www.codgenome.no/) Salmo salar data by using 
translated tblastn search (http://www.codgenome.no/blast/ 
blast_new.php) and the Trace Archive Nucleotide (http:// 
blast.ncbi.nlm.nih.gov/Blast.cgi). CCK-1R sequence frag- 
ments were identified in scaffolds gnl|ti|2263327594 and 
gnl|ti|2287331943 while CCK-2Rl/gastrin and CCK-2R2 
in scaffolds gnl|ti|2279893979, gnl|ti|2263 1015 19, and gnl| 
ti|2319219166gnl|ti|2317620741, respectively. ORF and 
5' untranslated regions were obtained by real-time poly- 
merase chain reaction (RT-PCR) and rapid amplification 
of cDNA ends (RACE) method according to the manu- 
facturer's instructions (Marathon™ cDNA Amplification 
Kit; Clontech, Mountain View, CA). Finally, cDNA frag- 
ments containing the complete coding regions of CCK-1R 
and CCK-2Rs were amplified using LongRange PCR Kit 
(Qiagen, Hilden, Germany). 

To amplify fragments within the range of 1.5 kb, PCR 
reactions were performed using GoTaq PCR reagents 
(Promega, Madison, WI), according to manufacturer's 
instructions. PCR conditions were as follows: initial acti- 
vation of the Taq-DNA-Polymerase for 5 min at 94°C, 
followed by 33 cycles of denaturation for 25 sec at 94°C, 
annealing for 30 sec at 58°C, and extension for 90 sec at 
72°C. The program ended with an incubation of 7 min at 
72°C. Gel-purified products were cloned into a pCR4- 
TOPO vector using TOPO TA Cloning® Kit (Invitrogen). 
Sequencing was performed at University of Bergen 
Sequencing facility using BigDye Terminator v3.1 chemis- 
try in ABI PRISM377 DNA sequencer (Applied Biosys- 
tems, Foster City, CA). 

Sequences were analyzed by homology across species 
with the basic local alignment search tool (blast; www.ncbi. 
nlm.nih.gov) to confirm gene orthology. Subsequent 
sequence analysis using other bioinformatics tools such as 
clustalW alignment revealed the finding of two paralogue 
duplicates (CCK-2Rl/gastrin and CCK-2R2) of the CCK- 
2R gene and one single copy of the CCK-1R gene. The lists 
of primers used to identify these genes are given in Table 1. 

Real-time quantitative PCR 

mRNA expression levels of CCK-1R, CCK-2R1 and CCK- 
2R2 and of an endogenous housekeeping gene encoding 
for Elongation factor 1 -alpha (Elf la; Genbank accession 
No. AF321836) were quantified using quantitative PCR 
(qPCR) analysis on the CFX-96 RT-PCR detection system 
platform (Bio-Rad, Hercules, CA), by Power SYBR Green 
PCR kit (Applied Biosciences, UK), in a final volume of 
25 /.iL per reaction. qPCR analysis was conducted on 
cDNA that was reverse transcribed using RNAs pooled from 
six fish. To minimize sample variability, all samples were 
run in triplicates. The qPCR conditions were as follows: 
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Primers 


Sequences 


Purpose 


cckar7 Fw 


GTACCTGGCAAACCAAGTCCCAT 


cckar first PCR 


cckar 9 Rv 


GATGGAGAAAGGAAGCTAAGCCTG 


cckar first PCR 


cckbr 5 Fw 


TCCACCTTCAGCCTGGTTGCCATA 


cckbr first PCR 


cckbr 9 Rv 


AGTCAGTCAATGGTTCCTCCCTCA 


cckbr first PCR 


cckar 10 Fw 


AGTGATGTCATCCAGCAGTCCT 


cckar RT-PCR 


cckar 1 1 Rv 


TCAGCTTGGACTTGTTGTTCTCTCT 


cckar RT-PCR 


cckbr 1 0 Fw 


TCCACCTTCAGCCTGGTTGCCATA 


cckbr RT-PCR 


cckbr 1 1 Rv 


ACGTAGCATCCATCACCGTCGT 


cckbr RT-PCR 


cckar 10 Fw 


AGTGATGTCATCCAGCAGTCCT 


cckar Q-PCR 


cckar 1 5 Rv 


CTGAAGGTAGCAGCCGTCGTTGT 


cckar Q-PCR 


cckbr ! 16 Fw 


ATCTCCACCTTCAGCCTGGTTGCCATA 


cckbr- 1 Q-PCR 


cckbr , 18 Rv 


ACACATGCGTGCGATGGTGATGTTGT 


cckbr- 1 Q-PCR 


cckbr 2 21 Fw 


TCATCGACGGCAACCACGTCCAAATT 


cckbr-2 Q-PCR 


cckbr 2 23 Rv 


AGGAAACCAGCACAGGAAGAACAGAGCA 


cckbr-2 Q-PCR 


Elf 1a Fw 


GAGAACCATTGAGAAGTTCGAGAAG 


Q- PCR reference 


Elf 1a Fw 


GCACCCAGGCATACTTGAAAG 


Q- PCR reference 


B'RACE ccka P, 


CATGGGACTTGGTTTGCCAGGTACG 


cckar B'RACE 


5'RACE ccka P 2 


GATTGAAGGTAGAGACACTCACCGAGAT 


cckar B'RACE 


B'RACE cckb P, 


TATGGCAACCAGGCTGAAGGTGGA 


cckbr B'RACE 


B'RACE cckb P 2 


AGAGCCCTATGGGAACGCATGGAA 


cckbr B'RACE 


cckbr Lr fw 


AGATCCACTCAACTCATGAACCCATACTG 


Long range PCR 


cckbr Lr Rv 


AGTCAGTCAATGGTTCCTCCCTCA 


Long range PCR 


cckar-2 Fw 


AAGAGAGAACAGCAAAAATGGAACCA 


Long range PCR 


cckar-2 Rv 


GATGGAGAAAGGAAGCTAAGCCTG 


Long Range PCR 


Cckbr2d fw 


GCACATAGCACGAATATAATGTGC 


Long range PCR 


Cckbr2f rv 


CCAGTCAGTCAATGGTTCCTC 


Long range PCR 



94°C for 5 min and 42 cycles at 94°C for 30 sec, 60°C for 
30 sec, and 72°C for 30 sec. Absence of primer dimers 
and nonspecific products was verified in every qPCR assay 
by melting curve analysis (temperature reading every 0.2°C 
from 60°C until 95°C). 

The choice of Elf la as a reference gene was based on a 
previous study, where amplification of this gene showed a 
steady-state level of expression among different segments 
of salmon brain and intestine (Murashita et al. 2009). 
Expression data of Elf la were used for normalization. 
Standard curves (twofold dilution series of cDNA) were 
generated for target and Elf la genes by plotting the cycle 
threshold (Ct) obtained in qPCR analysis versus the loga- 
rithm of input quantity of RNA and performing a linear 
regression (Bustin 2000). The threshold was consistently 
set for 0.10 and analyzed using CFX manager software. 
The data were exported in to Microsoft Excel for further 
analysis. Ct values of triplicates were processed using 
Q-gene software (Simon 2003). 

Phylogeny 

Phylogenetic analysis was performed by PhyML using 
software Phylogeny.fr (Dereeper et al. 2008) run with a 
maximum likelihood method (aLRT) statistical test 



(Jagerschmidt et al. 1998) of branch support to build a 
tree with 100 bootstrap replicates, based on available 
CCK-1R and CCK-2R-like sequences retrieved from 
Genbank (NCBI). Complete amino acid sequences for 
each protein were used to construct the phylogenetic tree. 
CCK-R-like sequences from invertebrates [Honey bee 
(Apis mellifera), Beetle (Tribolium castaneum), and Nema- 
tode (C. elegans)] were included as outgroups. 

In silico analysis of protein sequence, structure 
modeling, and ligand docking 

Putative transmembrane domains were predicted using 
TMHMM 2.0 server (Denmark), which is part of the 
Simple Modular Architecture Research Tool. Transmem- 
brane protein 2D topology was designed using RbDe tool 
of Institute of computational biomedicine, Cornell 
University (Skrabanek et al. 2003). 

Three-dimensional models for CCK-1R and CCK-2R1/ 
gastrin were created and viewed using Rasmol. Because 
the alignment between CCK-2Rl/gastrin and CCK-2R2 
gave 90% similarity including the position of functional 
motifs, we have only modeled CCK-2Rl/gastrin. PDB-ID 
2KS9 chain A determined using nuclear magnetic reso- 
nance was found to be a good template (Identity -25% & 
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Table 2. Bioinformatics tools and databases used for sequence and structural analysis. 



Purpose 



Programmes used 



Primary sequence analysis 



Putative transmembrane prediction 
Residual model 
3D Modeling 

Protein 3D visualization 
RMSD calculation 
Cysteine bridge analysis 
Protein motif analysis 
Ligand docking 

Interacting residues for ligand binding 
Phylogenetic tree building 



BLAST (www.ncbi.nih.nlm.gov) 

EMBOSS (http://emboss.sourceforge.net/) 

CLUSTALW (www.ebi.ac.uk/services) 

l-TASSER (http://zhanglab.ccmb.med.umich.edu/l-TASSER/) 
(http://www.cbs.dtu.dk/services/TMHMM/) 
RbDe (http://icb.med.cornell.edu/services/rbde/diagrams) 
Modeller (http://www.salilab.org/modeller/) 

SWISS-MODEL (http://swissmodel.expasy.org) 
Rasmol (http://openrasmol.org/) 
DALI (http://ekhidna.biocenter.helsinki.fi/dali_lite/start) 
WHAT-IF (http://swift.cmbi.ru.nl/servers/html/listcys.html) 
Regex engine inbuilt in PERL (http://www.perl.org/). 
HADDOCK (http://haddock.science.uu.nl/services/HADDOCK/haddock.php) 
LigPlot+ (http://www.ebi.ac.uk/thornton-srv/software/LigPlus/) 
PhyML http://www.phylogeny.fr/ 



Similarity -40%) for modeling using BLAST against pro- 
tein data bank (PDB). This template satisfied the addi- 
tional constraint of possible disulfide bond between the 
cystine residues in the studied protein sequence. The 
modeled structure with I-TASSER for CCK-1R and CCK- 
2Rl/gastrin has TM-score -0.56 ± 0.15 and -0.58 ± 0.14 
with RMSD score -9.9 ± 4.6 and -9.6 ± 4.6 A, respec- 
tively. SWISS-MODEL reported Anolea and Gromos 
scores in favorable negative range for most of the mod- 
eled residues with final energy lower than 3500 kj/mol in 
average. Furthermore, the Modeler showed mean DOPE 
score of -7200 and GA341 -0.55 on average. WHAT-IF 
analysis showed the possibility of cystine bridges between 
residue positions 119, 208 and 130, 213 for CCK-1R and 
CCK-2Rl/gastrin, respectively. Docking with natural ligand 
cholecystokinin-8 (PDB ID 1 D 6G chains B-CCK-8) 
was performed using HADDOCK and based on the 
HADDOCK scores, potential ligand interacting residue seg- 
ments were considered. For both CCK-1R and CCK-2R1/ 
gastrin the residues of extracellular loop 2 (ECL2) and 
ECL3 were found as the best candidates. Potential ligand 
interaction residues were predicted using LigPlot+. The 
sequences were aligned to find out the conserved potential 
residues and domains in Atlantic salmon. The software 
programs used for all the analyses in this study are listed in 
Table 2. 

Results 

Cloning and sequence analysis 

This study reports the identification of both subtypes of 
cholecystokinin receptors (CCK-1R and CCK-2R). Fur- 
thermore, we found duplicates of CCK-2R gene, named 



CCK-2Rl/gastrin and CCK-2R2. The accession numbers 
assigned by NCBI for CCK-1R is JX017294 and that for 
CCK-2Rl/gastrin and CCK-2R2 are JX017295 and 
JXO 17296, respectively. 

Using bioinformatics and RACE PCR approaches, we 
obtained full-length products for all genes studied here. 
The isolated product for CCK-1R cloned from pancreas 
cDNA, has an ORF of 1418 bases including a CDS 
encoding a translated protein of 461 amino acids with a 
molecular mass of 51.29 kDa. The ORF of CCK-2Rl/gas- 
trin cloned from midgut cDNA consisted of 1608 bases 
and comprised a CDS of 1359 bases encoding for a 
putative protein of 453 amino acids. CCK-2R2 tran- 
scripts cloned from brain tissues contained an ORF of 
1426 bases and encoded for a protein of 449 amino 
acids. Molecular mass of both CCK-2R receptors was 
predicted to be of -50 kDa. Sequence homology blast 
search confirmed gene orthology with CCK receptor 
genes of other species. 

Alignment of salmon CCK-1R showed 73% amino acid 
sequence similarity with zebrafish and Nile tilapia CCK- 
1R genes and 63% similarity with the human ortholog 
(Table 3). CCK-2R1 /gastrin and CCK-2R2 showed 90% 
amino acid sequence similarity among themselves and 
shared 70% similarity with zebrafish and Nile tilapia 
CCK-2R gene while about 50% with human and mouse 
orthologs (Table 4). 

Phylogeny 

The phylogenetic tree indicated that invertebrate CCK-R 
genes clearly form a separate clade with 100% boot strap 
value (clade C in Fig. 1). Both receptor types cluster together 
with the respective CCK subtype receptor orthologs in two 
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Table 3. Amino acid sequence identity (%) between newly identi- 
fied salmon CCK-1R and known CCK-1R orthologs. 



Salmon 
CCK-1 R 



Zebrafish CCK- 1 R 


73 


Nile tilapia CCK-1 R 


73 


Pufferfish CCK-1 R like 


80 


Frog CCK-1 R 


59 


Chicken CCK-1 R 


63 


Human CCK-1 R 


63 


Mouse CCK-1 R 


55 



Table 4. Amino acid sequence identity (%) between newly identi- 
fied salmon CCK-2R and known CCK-2R orthologs. 





Salmon 


Salmon 




CCK- 2R1/gastrin 


CCK-2R2 


Salmon CCK-2R1/gastrin 




90 


Salmon CCK-2R2 


90 




Zebrafish CCK-2R 


70 


68 


Nile tilapia CCK-2R 


71 


71 


Pufferfish- CCK-2R like 


70 


70 


Frog CCK-2R 


55 


55 


Chicken CCK- 2R 


55 


56 


Human CCK- 2R 


49 


50 


Mouse CCK-2R 


47 


47 



main clades with good boot strap support (clades A 
and B in Fig. 1; 72 and 98%, respectively). Interestingly, 
CCK-2Rl/gastrin and CCK-2R2 genes are more closely 
related to each other (boot strap support 100%) than 
either of them is to their orthologs as predicted if they 
arose from recent lineage-specific genome duplication. 
These results are in good agreement with the concept of 
traditional taxonomy. 

Structural analysis 

Residual modeling of salmon CCK receptors clearly sug- 
gests that these proteins belong to seven-transmembrane 
GPCR family (Fig. 2). The 3D models generated taking 
PDB ID 2KS9 as template confirmed this assignment 
(Fig. 3). WHAT-IF analysis showed the possibility of 
disulfide bond among cystine residues between first and 
second extracellular loop (Fig. 4), further confirmed with 
prediction from ScanProsite database. The amino acid 
sequences were aligned with zebrafish, Nile tilapia, rat, 
mouse, and human to identify potentially conserved 
functional motifs. Functional motifs such as E/DRY, 
NPxxY, PLC, PKC, CRAC, CCM, and Palmitoylation 
domains (Kristiansen 2004; Potter et al. 2012) were 
found conserved in residual positions to the human 
CCK receptors (Fig. 5A and B). A detailed account of 
motif location corresponding to human and its function 
is given in Table 5. Sequence alignment confirmed the 
position of human CCK-8 affinity determinant residues 



-c 



-HoneybeeCCK-RJikeXP-003250082 
-Beetle_CCK-R_like_XP-972750 
-Nematode_CCK_BR_like_NP_00 1 022844 
-SalmonCCK-ARJXO 1 7294 
-Tilapia CCK-AR XP-00344942 1 
-Pufferfish_CCK-AR_like_CAG04576 
-ZebrafishCCK-ARXP-697493 
-Frog_CCK-AR_XP-002939147 
-Chicken_CCK-AR_NP-00 1 074970 
-Turkey_CCK-AR_XP-003205928 
-Lizard_CCK-AR_XP-003226175 
-Mouse JXK-ARJflP-033957 
-Monkey_CCK-AR_EHH25794 
-Human_CCK-AR_NP-00072 1 
-Salmon CCK-BR1 JX017295 
-Salmon_CCK-BR2 JX0 1 7296 
-Tilapia_CCK-BR_XP-003452984 
-PufferfishCCK-BR Jike_C AG05 857 
-ZebrafishCCK-BRXP-00266336 1 
-Frog_CCK-BR_NP-00 1 079277 
-Chicken CCK-BR NP-00 1 00 1 742 
-Mouse_CCK-BR_NP-03 1 653 
-Monkey_CCK-BR_EHH23 146 
-Human CCK-BR NP-795344 



Figure 1. Phylogenetic analysis of CCK-1 R (CCK-AR) and CCK-2R (CCK-BR). Numbers at nodes indicates the bootstrap value, as percentages, 
obtained for 100 repetitions. Gene accession number for each species is given along with the species name. Branch A, B and C represents 
vertebrate CCK-1 R, vertebrate CCK-2R, and invertebrate CCK-R-like genes, respectively. 
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Figure 2. Residual model diagram of salmon CCK-1R (A) and CCK-2R1 (B). 



Leu- 116, Phe-120 of first ECL; Tyr-61of TMD1 (trans- 
membrane domain 1); Met 186, Tyr 189, Thr 193 of 
TMD3; His 207 of second extracellular loop were con- 
served for both the salmon receptors studied here 
(Fig. 5). Ligplot+ docking analysis data further sup- 
ported possible interaction of unsulfated CCK-8 with 
Phe-125 in first extracellular loop and His-377 in TMD7 
for CCK-2Rl/gastrin (Fig. 6). 

Gene expression and tissue distribution 

The mRNA tissue distribution of the three CCK receptors 
in eight tissues from Atlantic salmon was assessed using 
qPCR approach (Fig. 7). CCK-1R mRNA expression was 
significantly abundant in pancreas followed by midgut, 
hindgut, and gallbladder. Relatively low levels of CCK-1R 
transcripts were found in pyloric sphincter and pyloric 
ceca while the lowest levels were found in brain. A signifi- 
cantly higher level of CCK-2Rl/gastrin was seen in 
midgut compared to the other organs. The levels of CCK- 
2Rl/gastrin in stomach, gallbladder, and pancreas were 
not significantly different from each other. CCK-2Rl/gas- 
trin mRNAs were not detected in pyloric sphincter and 
hindgut. CCK-2R2 mRNA expression was most abundant 
in brain followed by exocrine pancreas, pyloric sphincter, 
and hindgut. CCK-2R2 mRNA transcripts were not found 
in stomach and gallbladder. 



Discussion 

Sequence and phylogeny 

The predicted evolutionary relationship of CCK-R and 
CCK-R-like receptors was evaluated using phylogenetic 
analysis with annotated ortholog sequences from the 
Genbank database. In our phylogenetic tree, vertebrate 
and non vertebrate CCK-R fall into separated clusters 
and a clear evolutionary split between the CCK-1R and 
CCK-2R is observed in the vertebrate group (Fig. 1). 
This is consistent with Salmon CCK-2Rl/gastrin and 
CCK-2R2 showing 47 and 49% amino acid sequence 
homology with CCK-1R, respectively. Thus, although 
CCK-1R and CCK-2R genes are predominantly con- 
served among species, they show important functional 
divergences and exhibit different affinities for gastrin and 
CCK ligands (Yassin 1999). Besides, positions of both 
receptors in the tree pinpoint Nile tilapia as the closest 
relative to salmon. However, because only a limited 
number of teleost CCK-R sequences are currently avail- 
able, this aspect remains to be explored further. 

Duplicated genes are frequently found in teleosts as 
result of the fish-specific whole-genome duplication 
(FSGD, 3R duplication) that occurred about 350 million 
years ago (Meyer and Schartl 1999; laillon et al. 2004). 
An additional genome duplication event was postulated 



2013 | Vol. 1 | Iss. 5 | e00069 
Page 6 



© 2013 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of 
the American Physiological Society and The Physiological Society. 



R. M. Rathore et al. 



CCK Receptors in Atlantic Salmon 





Figure 3. Rasmol visualization of 3D model of Atlantic salmon CCK-1R (3.1) and CCK-2R1/gastrin (3.2) Each figure is shown in four panels: 
(a). Top view (b). Front view (c). Front view with predicted motif positions (d). CCK-8 and CCK-R complex. 
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(A) 

1 48 

MEPFTLHDMLINSTKLYKILCDFGIKHVSECEDERETPPEPKDLNQTV RIFLYSLIFLVSVLGHSLIIAVLV RNRRMRTVTNL 

TMDl 

92 108 112 119 143 160 

FLLSLAASDLMLCLFCMPFTLIP NLMRDFVFGSGI©C VAMYFMGI5V5VSTFNLVAI 5l )ERY^ AICHPLTSRTWQTK SHAAKV 
TMD2 j TMD3 

2K)8 236 242 246 

ISATWWSFLFMLPYPI SS TLVPFTRVMMSTGNI^tLVWPSDVIQQSW YVSLLLLLFLVPGIMMMTAYGLI SLELYRG I KFEM 
TMD4 TMD5 

325 

ANRKSSRERQCSTGSIKPGDNDGCYLQPAKKKGELAHLQNPPSTPSPNMTINSRESNKSKLSRVCSNSSTCNLMA^KRj/IRML 

340 350 383 404 

LVIVCLFFLCWTPVFAVMAWRA FDRRSADKLLSGAP ISFIHLLSYTSAC^MPIIY^FMH KRFRQGILSTFTfcC^SPKAGGVGR 
TMD6 TMD7 

451 

GAGNRMGTGRGGGRGGMRSGEENGHMPQSGMl f ^ CA 

(B) 

60 64 

MDLQTLMETSVFLECVEAQGAASHGTPNSSLELTNITTCGNAYWLKPQHVKQKEVDSLRILLy SVIFLLSVFGHLLIIWLT 
MDLRTLNETSTFLECVGAQGAASNGTMNVSLGLTNITTCGNESWLKPQHVKQKEVDSLRILLYSVIFLLSVLGSLLIIWLa 

TMDl 

95 103 125 130 154 

VMKRMRTVTHSFLLSLAVSDLMMA IFCMPFTLIP NLLEDFIFGPAMgyi VAYLMGVSVSISTFSLVAIAllERVlp AICHPLKSR 
VBKRMRTVTNSFLLSLAVSDLMMAIFCMPFTLIPNLLEDFIFGPAMCKIVAYLMGVSVSISTFSLVAIAIERYSAICNPLKSR 
TMD2 I TMD3 

171 178 189 196 213 235 254 

AWQTRSH AYRVITATWLLSFMIMSPYPVFS HLVHVPLKDNITIARMg^HIWPHREVEQTWN MMLLLTLFWPGWMIVAYGLI 
AWQTRSHAYRVIAATWLLSFMIMTPYPVFSHLKHVPLKDNMTITRMCRHIWPYPEVEQTWNMMLLLTLFWPGWMIVAYGLI 
TMD4 TMD5 

267 329 

srelyrgiqfelgqktsspglkngltgtvscgsddgdgcyvqvskrphsmemstltsstastskveharsntseakll/{kk^v 
srelyrgiqfemvqktsspglkngltttgacgsddgdgcyvqlakrprsmemstltsstattsklerprsntseakllakkrv 

343 347 354 377 387 409 

IRMLVVIVAL: FLC MPLYCANTWKAFHPASARRALSG APISFI H LLSYTSAC\|NP~ I Y^ FMNTRFRQSLLATFACjCCfcPPCR 
IRMLWIVALFFLCWFPLYCAKTWKAFHPVSARRALSGAPISFIHLLSYTSACV1]?ITYCFMMTRFRTALLATFACCCRPPCR 
TMD6 THD7 
437 

H — QGLRDGEEDAMALGVSM jKFSYTTVST^ flGRAEGGTID 
HRRHGLRDGEEDAMAMGVSMSKFSYTTVSTMGTC 

Figure 4. Predicted functional motifs and affinity determinants in the amino acid sequence of Atlantic salmon CCK-1R (A), and CCK-2R1/ 
gastrin and CCK2-R2 (B) CCK-2R1 sequence is presented on the top and CCK-2R2 below. Functional motifs are marked red with each residual 
position given above the sequence. Trans-membrane domains are underlined and Cystine bridge position is circled (for description see Table 5). 



to occur in the last common ancestor of salmonids (4R 
duplication), resulting in a pseudotetraploid karyotype 
with multiple gene copies (Ohno 1970; Moghadam et al. 
2005; Kurokawa and Murashita 2009). CCK-2Rl/gastrin 
and CCK-2R2 paralogues share rather high amino acid 
conservation (90% similarity) suggesting that they most 
probably are recently duplicated genes following the 4R 
duplication of salmonids (Table 4). Furthermore, we 
postulate that salmon CCK-2Rl/gastrin is the ortholog of 
mammalian gastrin receptor due to its prominent expres- 
sion in alimentary-tract tissues and also due to the fact 
that CCK-2R2 is clearly predominant in brain. 

Structure and motif analysis 

CCK receptors belong to the /?-group of rhodopsin recep- 
tors or to the subfamily A6 in the A-F classification of 
GPCRs (Staljanssens et al. 2011). The residual and 3D 
structural modeling of salmon CCK-1R and CCK-2R1/ 
gastrin was predicted to be typical of the GPCR superfam- 
ily by having seven TMDs (Fig. 3). We have verified the 
conventional characteristic signature sequences in human 
that provide the appropriate architectural motifs, to 
achieve the preferred structure for proper activity and reg- 
ulation of the receptor. The alignment of CCK-1R and 



CCK-2Rl/gastrin with respective ortholog sequences 
(Fig. 5A and B) clearly suggest that all these functional 
motifs, including their respective positions, are conserved 
in salmon for both the receptors. This is also the case for 
the signal transduction motifs PLC and PKC (refer to 
Table 5 for details). 

Ligand affinity 

We were also interested to analyze ligand interacting 
regions in the sequences of salmon CCK-R genes. We 
have considered important residues and their positions 
for ligand affinity based on human and cell based studies. 
It is well known that CCK receptors provide a common 
structural basis for affinity to their natural ligand CCK, 
known to be present in several molecular variants. In 
humans the octapeptide CCK-8 (Asp-Tyr(S0 3 H)-Met- 
Gly-Trp-Met-Asp-Phe-NH 2 ) has the highest activity (Escri- 
eut et al. 2002). CCK-1R has a greater affinity for sulfated 
CCK than for sulfated or unsulfated gastrin and requires at 
least the C-terminal heptapeptide as requisite for binding 
and to elicit biological activities. On the other hand, CCK- 
2R/gastrin also binds to gastrin since this has the same C-ter- 
minal tetrapeptide and does not distinguish between gastrin 
and CCK, whether sulfated or not (Shulkes and Baldwin 
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(A) 



(B) 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 

Rat CCK-AR 
Mouse CCK-AR 



Salmon CCK-AR 
Zebrafish CCK-AR 
Tilapia CCK-AR 
Hat CCK-AR 
Mouse CCK-AR 
Human CCK-AR 



Salmon CCK-AR 
Zebrafiah CCK-AI 
Tilapia CCK-AR 
Rat CCK-AR 
Mouse CCK-AR 
Human CCK-AR 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 

Rat CCK-AR 
Mouse CCK-AR 



:-ar 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 
Hat CCK-AR 
Mouse CCK-AR 
Human CCK-AR 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 
Rat CCK-AR 
Mouse CCK-AR 
Human CCK-AR 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 



Salmon CCK-AR 
Zebrafiah CCK-AR 
Tilapia CCK-AR 

Rat CCK-AR 
Mouse CCK- AH 
Human CCK-AR 



Salmon CCK-AR 
Zebrafish CCK-AR 
Tilapia CCK-AR 



MEPFTLHDMLINSTNLYKILCDFGIKNVSECEDERE-TPPEPK-- 42 

METFTIQDMLVNSTDIYKILCGLGIGNLSECKNESTPSPSEPK-- 4 3 

METFTIHDMLVNSTDLNRILCSFGIKNISECESEQD-LL PEPKAY 44 

MSHSPARQHLVESSRMD- - WDSLLMNGSNIT- PPCELGLENETLFCLDQP -QPSK 52 

MD--WDSLLMNGSNIT-PPCELGLENETLFCLDQP-QPSK 3 7 

MD--VVDSLLVNGSNIT-PPCELGLENETLFCLDQP-RPSK 3 7 

-- DLNQTVRIFLYSLIFLVSVLGNSLIIAVLVRNRR 76 

- - DMNQTVRIVLYALIFILSFLGNSLI IAVLLRNRR 7 7 

GLCIHLNKCRNWYLCACKLVLCPTSTDINQTVRIALYSLIFLLSVLGWSLIIAVLVRNRR 104 

EWQSALQILLYS I IFLLSVLGNTLVITVLIRNKR 86 

EWQSAVQILLYSF IFLLSVLGNTLVITVLIRNKR 71 

-- - - EWQPAVQILLYSL IFLLSVLGNTLVITVLIRNKR 71 



MRTVTNLF LL S LAAS D LMLC LFCM PFTLI PMLMR D F VFG SG Tjj^KVAMYFMG I SVSVSTFN 136 

MRTVTNLFLLSLAVSDLMLCVFCMPFTLIPNLMKDFIFGSGM&CVATYFMGISVSVSTLN 137 

MRTVTNLF LL S L AVSD LMV S LV C I PFTL I PNLMR DF I FGNG T@CL VMY FMG VS V S I ST FN 164 

MRTVTNIFLLSLAVSDLMLCLFCMPFNLIPNLLKDFI FGS A\©CTTT Y FMGTSVSVSTFN 14 6 

MRTVTNIFLLSLAVSDLMLCLFCMPFNLIPNLLKDFIFGSAV@CTTTYFMGTSVSVSTFN 131 

MRTVTNIFLLSLAVSDLMLCLFCMPFNLIPNLLKDFI FGS AV@CTT T Y FMG TS V S VS T FN 131 



S LE R YS A I CNPL TS RTHQ TK SHAAKV I S AT WWS F LFM LPYPISSTLVPF TR VNNS 196 
SLEHYSAICNPLKSRAWQTKSHAVKVISATWWSFLLMLPYPIGSTLVPFIRSNNS 197 
SLERYSAI CNPLTSRAWQTKSHAAKVI TATWLGS F I LMLPYP I FSTLKPFTRLNNS 224 
SLERYGAICRPLQS R VWQTKS HAL KV I AAT WC L S FT I M T P Y P I Y SNL V P F TKNNNQ 2 06 
SLERYGAICRPLQS R VWQTKS HALKVIAATWCLSFTIMTPYPIYSNLVPF TKNNNQ 191 
SLEHYGAICKPLQSRVWQTKSHALKVIAATWCLSFTIMTPYPIYSNLVPFTKNNNQ 191 



tgnmsrlwpsdviqqswH 



'LVPGIMMMTAYGLI SLELYRGIKFE MANR KS SR 2 56 
TGNMgRLVWPSDVMQT-WYIFLLLTLFLVPGIIMMTAYGLTSLELYRGIRFEMTNRKSNR 2 56 
TGH»©iLVWPNDIIQQSWYVSLLLLLFLIPGIVMMTAYGLISIELYRGIKFEMSKRKVGR 2S4 
TANMjSiFL L P SDAMQQ SWQTFLLLILFLLPGI VH W AY GLISLELYQGI KFDASQKKS AK 266 
T ANM©tFL L P SDAMQQ SWQTFLLLILFLI PGWMWAY GLISLELYQGI KFDASQKKS AK 251 
T ANM§JiFL L PNDVMQQ SHHTFLLLILFLIPGI VMMV AY GLISLELYQGIKFEAS QKKS AK 2 51 



ERQCSTG- - 



-sikpgdndgcylqpakkkgelahlqnppstpspnmtinsee;:- , :;k ■ ■ - 



EKE Y GTS SLKNGDSDGCYFQPSKKK- - 

DRQSSTS SIKPGDSDGCYLQPSKKK- - 

EKKPSTG SSTRYEDS DGC Y LQKSRPPR- 

EKRL S SGGGGGGG SSSSRYEDSDGCYLQKSRPPR - 
ERKPSTT SSGKYEDS DGC Y LQ KTRP PR - 



-RSEPVMPS HPKS 293 

SLTMNT GSS 31B 

-KLELQQLS SGSGGS 307 

-KLELQQLS TSSSGG 29S 

■ - KLELRQLS TGSS-S 291 



N -SSTCNLMAKKRVI RMLLVIVCLFFLCWTPVFAVNP 

KLRRVCSN-SPTANLMAKKRVIRMLLVIVGLFFLCWTPVFWNAWRAFDRRSADSLLSGA 3 52 
MVGR VCGSGS ST SNLM AKKRV I RMLL V IVFLFFLCWTPV FWNA WQ AFDRR S A Y R - LTGA 3 77 
RLNR IRSS - S SAANLIAKKRVI RML I VI WLFFLCWMP I FS ANAWRAYDTVS AEKHLSGT 356 
RINRIRSS-GSAANLIAKKRVIRMLIVIWLFFLCWMPIFSANAWRAYDTVSAEKHLSGT 358 
RANR I RSN - S SAANLMAKKRVI RMLI VI WLFFLCWMP I FSANAWRAYDTASAERRLSGT 350 



PISFIHLLSYTSACVNPIIYCFMNKRFRQGILSTFTCCSSPKAGGVG HGAGNRMG A 

PISFIHLLSYTSACVNPI VYCFMNKRFRQGILSTFSCCRG-DVDEFS RSSTRRSA ' 

PISFIHLLSYTS ACVNP 1 1 YCFMNKRFRQGMLAT FSCCSCFR TGGGG S S S L TR SAG S R LT A 

PISFILLLSYTS SCVNP I I YCFMNKRFRLG FMAT F PCC PN PGPPGV A 

PI SFI LLLSYTS SCVNP 1 1 YCFMNKRFRLGFMATFPCCPN - - PGPTGV A 

PISFILLLSYTS SCVNP IIYCFMNKRFRLGFMATFPCCPN PGPPGA : 



TGRGGGRGGMR- SGE ENGHM PQ SGNNT RFTYSSIRGSAQA- - - A 

RGAGGLGT LFKGAGQTDGHTQT SGTLTKLTDTSVRGSAQA • 

KGEAGRPREGPKSPEQNGHTSPSGTSTRFTYTGIRASAWSELT A 

RGEVG EEEDGRTIR-ALLSRYSYSHMSTSAPPP- - '. 

RGEVG EEEDGRTIR-ASLSRYSYSHMSTSAPPP - - A 

RGEVG EEEEGGTTG-ASLSRFSYSHMSASVPPQ-- A 



Salmon CCK-BR1 
Salmon CCK-BR2 

Tilapia CCK-BR 
Zebrafiah CCK-BR 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 

Tilapia CCK-BR 
Zebrafish CCK-BR 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 

Zebrafish CCK-BR 

Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 

Tilapia CCK-BR 
Zebrafish CCK-Bi 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 

Salmon CCK-BR2 
Tilapia CCK-BR 
Zebrafish CCK-BI 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 
Tilapia CCK-BR 
Zebrafish CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 
Tilapia CCK-BR 
Zebrafish CCK-B! 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



Salmon CCK-BR1 
Salmon CCK-BR2 
Tilapia CCK-BR 
Zebrafish CCK-BI 
Rat CCK-BR 
Mouse CCK-BR 
Human CCK-BR 



MDLQTLNETSVFLECVEAQGAASNGTPNSSLELTNITTCGHAYWLKPQHVKQKEVDS - L 59 
MD L RT LN E TS T FLE C VGAQGAASNGTMNVS L GLTN I TTCGNE S WLK PQHV KQKE VDS - L 59 

MPPFYTDEPSS SYCEDCPEG TDEDHT-L 27 

MDLQLLNDTGIWACLQLHNASILKRSSNGSLE-VPNGTCEDALLLPKIEPARAKEMDS-L 58 

MEL LKLNRS VQG PGPGSGS S LCR PGVS LLNS SS AGNLS C E P P R I RGTGTRE LEMA I 56 

MD LLKLNRS LC.G PGPGSGS S LCR PGVS LLNS SS AGNLSC ET P RIRGTGTRELELTI 56 

MELLKLNRS VQGTGPGPGAS LCR PGAP LLNS SS VGNLS C E P P R IHGAGTRE LE LA I 56 



EILLYSVIFLLSVFGNLLIIVVLr.7"-:R.'-:FTVTn^ : LL: : LAVS" LMMAIFCMPFTLIPNL 11 S 

RILLYSVIFLLSVIGSLLI IVVLAVNKRMRTVTNSFLLSLAVSDLMMAIFCMPFTLIPNL 11 S 

RILLYSVIFFLSVFGNLLI IVVLTVNKRMRTVTNTFLLSLAVSDLMMAIFCMPFNLIPS I B7 

RILLYCLIFLLSVFGNLLI IWLMVNKRMRTVTNS FLLSLAISDLMMAWCMPFTLIPNL 118 

HITLYAVIFLMSVGGNVLI IWLGLSRRLRTVTNAFLLSLAVSDLLLAVACMPFTLLPNL 116 

RITLYAVIFLMSVGGNVLIIWLGLSRRLRTVTNAFLLSLAVSDLLLAVACMPFTLLPNL 116 

RITLYAVI FL?'[5V 7 7I7-.7 ~ 7.7'7 7777 77?.7\™7--F77SL.-.VSD77LAY.-i7'[''PF717FiJL II 'J 




IMG VSVS I STF S LVA IAIEF.YSAI CN P LKS R AWQT R S H AYRV I T 17 S 

IIVAYLHGVSVSISTFSLVAIAIERYSAICNPLKSRAWQTRSHAYRVIA 17 3 
IVTYLMGISVSISTFSLVAIAIERYSAICNPLKSRVWQTRSHAYRVIA 14 7 
I VA YLMG I SVC I STF S LVA I AI E RYS A I CN P L KS R AWQT R SHAY KV I A 17 3 
! S Y LMG VSV S VSTLNLVA I ALERYS A I CRPLQAR VWQT RSHAARVI L 17 6 
.VS Y LMG VSVS VSTLNLAA I ALE RYS A I CR P LQAR VWQT R SHAARV I L 17 6 
.VS YLMG VS VS VSTL S LV A I ALERYS A I CRPLQAR VWQT R SHAARV I V 176 



ATWLLSFMIMSPYPVFSH LVHVP LKDN I TI ARMi 
ATWL L S FM I MT P Y P V F S HL KHVP L KDNMT I T RMi 
ATWVLAL I IMI PYPI ISHLESFPRRDNTTAHR-< 
STWWSFLIMT PYPI FSTLVEFT-KNNNITGW 



IHIWPHREVEQTWt 

(H I W P Y PE VEQTWNMML LLTL FWP G 239 

! H KW P VAT AEQTW Y ILLLLVLFAIPG 206 

(HDWKNSEVEKAWCILLLFVLFFIPG 237 



ATWLL SG LLMV P Y P V YTMVQ PVG PRV LQ@1HRW P S AR VQQTW SVLLLLLLFFIPG 231 

ATWLL SG LLMV P YP V YTWQPVG PRILQ&1HLWPSERVQQMWSVLLLILLFFIPG 231 

ATWLL SG LLMVP YP VYTWQP VG PRV LQ@/HRWPSARVRQTW SVLLLLLLFFIPG 231 



WM IVAYGLISRELYRGI QFELGQKTS S PGLKNGLTGTVS C 2 8 C 

WM I VAYGL I S RELYRG IQ FEMVQKTS S PGLKNGLTTTGAC 28 C 

VMM I VAYGL I S RE LYRG IQ FEMGHKKD S SDVN - GLTATVSA 246 

WMIIAYGLISRELYRGIQFELEQKKGQ CGMKNG I S ATVL I 27 S 

W I AVAYGL I S RE L Y LGLH FDG E ND S ETQS RARNQGGL PGG - AAPG P VHQN GGCRPVTSV 29C 

WMAVAYGL I SRELYLGLR FDGDNDS ETQSRVRNQGGL PGGAAAPGP VHQNGGCRHVTS L 291 

WMAVAYGL I S RELYLGLR FDGDSD SD SQSRVRNQGGL PGA VHQNGRCR P E TG A 28 5 



GSDDGDGCWQVSKRPHSMEMSTLTSSTASTSKVEHARSNTSEAKLLAKKRHB ,I,VVIV 3 A - c 

GSDDGDGCYV Q LAKR PRSMEMSTLTS S TATTSKLE RPRSNT S EAKLL AKKRV I RMLW I V 34 C 

GSD DGDGCYVT AVQR PHTM EMS AMAAS - GRA I KAE RPRSNT S E AKLE AXKR V I RMLW I V 305 

SSDDGDGCYIQVSKRPNSMEMSALTPT--GSAKVDRPRSNTSEAKLMAKKRVIRMLIVIV 336 

AGEDSDGCCVQLPRSR- -LEMTTLTTPTPGPVPGPRP NQAKLL AKKH WRML LV I V 344 

TGE D SDGCYVQ L PRS R - - L EMTTLTT PTTG PG PGPRP NQAKLL AKKR WRM L LV I V 345 

VGEDSDGCYVQLPRSRPALELTALTAP--GPGSGSRP TQAKLL AKKR WRMLLV I V 33 9 



ALFFLCWMPLYCANTKKAFHPASAP.HALSGAPISFIHLLSYTSACVNP 1 1 YCFI'INTRFRQ 40C 

ALFFLCWFPLYCANTWKAFHPVSARRAL-3GAPISFIHLLS YTSACVT]? " 7 YCFJ4NTRFRT 40C 
VLFFVCWMPLYCVNTWRAFDESSATQTLSGAPIAFIHLLCYTSACVNPIIYCFMNTRFRK 365 
AMFF I CWMPVYS VNTWKAF DLAS AHRALSGA P I S F I HLL S YT S ACVNP 1 1 YC FMNKRFR K 39 6 
LLFFLCWLPVYS VNTWRAFDG PGAQRA LSGAP I S F I HLL S YVS ACVNP LVYCFMH RRFRQ 404 
LLP FVCW LPVYS ANTWRAFDG PG ARRA LAG A PISFIHLLSYT S AC ANP LVYC FMH RRFRQ 40 5 
VLF FL CW LPVYS ANTWRAFDG PGAHRALS G A P I S F I HLL S YAS ACVNP LVYCFMH RRFRQ 39 9 



■3li.^77a: .-:777 :7.^- -;;i?.:7-zi:a77-.l;v;:-:7 ■. ••' • ' •• 

al latfacccr p pcrhrrhglrdge ed amamgv smskfsyttv stmgtc 44 9 

al lat fs cc i apchrrrhrglrdne e d amat g asm s kfs yttv etmgtc 414 

al lat fs cccw pcrkr gfregdde vtatgasm s kfnyttvstmg pc 442 

ac ldtcarcc p r p prarpq plpd-edpptpsiasl srls ytt i stlg pg 4 52 

acldtcarccprpprarprplpd-edpptpsiaslsrlsytti stlgpg 453 

acletcarccprpprarpralpd- edpptps i aslsrlsytti stlgpg 447 



Figure 5. Alignment of CCK-1R (CCK-AR) amino acid sequences (A), CCK-2R1/gastrin and CCK-2R2 (CCK-BR) (B) with known functional 
domains of other vertebrates. Trans-membrane domains are marked in green. Residues marked bold represent functional motifs and those in 
circles cystine bridge residues. 



1997; Dufresne et al. 2006; Staljanssens et al. 2011). Cell 
mutational studies suggest that all the amino acids deter- 
mined to affect CCK-8 affinity in the CCK-2R protein 
are not necessarily important for CCK-8 binding in the 
CCK1R. However, two residue positions (Leu and Phe) of 
ECL1 seem to affect CCK-8 affinity in both CCK-2R and 
CCK-1R (Silvente-Poirot et al. 1998). These residues are 
conserved in salmon for both CCK-1R and-2R proteins 
(Fig. 5AandB). 

Human CCK-1R residues in ECL 2 (Met and Arg) and 
TMD6 (Arg and Asn) that account for the selectivity and 
receptor-ligand stability for sulfated versus nonsulfated 
ligand, are found conserved with salmon. Interestingly, 
other reported signature residues for sulfated analogue 
binding (Trp39, Gln40) in human CCK-1R that interact 
with N-terminal moiety of CCK (Kennedy et al. 1997) are 
not found in salmon and suggest possible differential 
N-terminal binding mechanisms for sulfated CCK in this 
species. This is also reflected in the Ligplot+ analysis for 



CCK-1R with unsulfated CCK-8, where the prediction 
plot did not show any significant ligand-receptor 
interaction as per agreement with the available literature 
(Fig. 8A and B). This suggests that salmon CCK-1R may 
also have a preference for sulfated analogs over unsulfated 
ones as that in mammals. However, these aspects are sub- 
jected to further conformational studies. 

For CCK-2R, experimentally determined putative CCK- 
8 affinity determinants in TMD 1 (Tyr), TMD4 (Met, 
Tyr, Thr), TMD 6 (Asn), and residues in ECL2 and ECL3 
are all conserved with the position of human CCK-2R 
(Silvente-Poirot et al. 1998; Anders et al. 1999; Gigoux 
et al. 1999; Blaker et al. 2000; Gouldson et al. 2000; 
Giragossian and Mierke 2001; Gale et al. 2003). The 
Ligplot+ analysis of salmon CCK-2Rl/gastrin predicts 
interaction of CCK-Tyr with salmon His-377 in TMD7 
(Fig. 6). This is in agreement with earlier reports based 
on site-directed mutagenesis data for 3D model of the 
human CCK-2R and CCK complex (Foucaud et al. 2008). 
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Figure 6. Ligplot+ predictions of first Extracellular loop (A) and second Extracellular loop (B) of Atlantic salmon CCK-2R1/gastrin. Figure shows 
complex of unsulfated CCK8 (molecular model in dark blue) with interacting residues (marked in green). 
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Figure 7. Tissue distribution and relative expression levels of the 
CCK receptors in Atlantic salmon. CCK receptor mRNA levels were 
normalized with Elf 1a gene. Error bars represent standard error of 
the mean (n = 3). Br, brain; St, stomach; Ps, pyloric sphincter; Gb, 
gall bladder; Pa, pancreas; Pc, pyloric ceca; Mg, midgut; Hg, 
hindgut. 

Another study on human models reported that CCK-Trp 
binds to Phe-120 located in the hydrophobic/aromatic 
pocket formed by residues from TMD3,7 and from the 
first extracellular loop (Foucaud et al. 2008). The align- 
ment of amino acid sequence suggests the conservation of 
residual position between human Phe-120 and salmon 
Phe-123. Our ligand docking analysis from Ligplot+ pre- 
dicts that Phe-123 in the first extracellular loop of salmon 



CCK-2Rl/gastrin interacts with CCK-Tyr while residue 
Phe-125 with CCK-Trp (Fig. 6). However, salmon Phe- 
125 is also found to be conserved in human Phe-122 
(Fig. 5B). As these data are based on predictions, further 
biological studies are essential to establish the exact inter- 
acting positions. 

Tissue localization and expression pattern 

The tissue expression pattern for the CCK receptor genes 
in Atlantic salmon differed substantially. We found that 
CCK-1R and CCK-2Rl/gastrin had a higher expression in 
tissues connected with direct regulation of gastrointestinal 
events than CCK-2R2 and that CCK-2R2 had its highest 
expression in brain. The spatial expression patterns of 
different salmon receptors analyzed in this study can sub- 
stantiate the hypothesis of functional conservation during 
evolution of these molecules with their ortholog counter- 
parts. Murashita et al. (2009) have reported the mRNA 
distribution of two CCK peptides in Atlantic salmon 
namely CCK-N and CCK-L. Phylogenetic analysis sug- 
gested that CCK-N and CCK-L belonged to the CCK-1 
and CCK-2 gene families respectively. Interestingly, 
salmon CCK receptors also follows similar evolutionary 
patterns. 

CCK-1R is reported to be involved in satiation via links 
to peripheral nervous system and vagus nerve. The most 
accepted vertebrate scheme is that the nutrients enter duo- 
denum and stimulate CCK secretion, some of the CCK acts 
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Figure 8. Ligplot+ prediction of 1 st Extracellular loop (A) and 2 nd Extracellular loop (B) of Atlantic salmon CCK-1R. Figure shows complex of 
unsulfated CCK8 (molecular model in dark blue) with interacting residues (marked in green). 



in a local paracrine manner to stimulate CCK-1R on the 
sensory fibers of the vagus nerve (Woods 2004). In mam- 
mals, CCK-1R is involved in both endocrine and exocrine 
functions of pancreas, in the regulation of GI motility and 
gallbladder contractions and it controls the release of 
enzymes such as somatostatin and pepsinogen from gastric 
mucosa (Lee et al. 1993; Crawley and Corwin 1994; Wank 
1995). Available data suggest that CCK-8 released in the 
CNS after a meal is responsible for the postprandial 
increase in colonic motility and that these effects are medi- 
ated through activation of central CCK-1R (Gue and 
Buenno 1991; Gouldson et al. 2000; Lin et al. 2000). A 
recent study on yellow-tail fish reported that the expression 
pattern of CCK-1R in gut is similar to that of humans 
(Furutani et al. 2013). Furthermore, these authors observed 
that CCK stimulates CCK-1R mediated gall bladder con- 
traction. In this study, the mRNA expression profile of 
CCK-1R gene in gut tissues of Atlantic salmon are consis- 
tent with these findings. 

CCK receptors act as gateway of diverse cellular event 
such as CNS-related functions, growth, proliferation, and 
inflammatory responses depending on their location and 
metabolic states (Nagata et al. 1996; Yassin 1999; Piiper 
et al. 2003; Luyer et al. 2005). In brain, CCK-2R has been 
demonstrated to participate in cognitive process, release 
of neurotransmitter and other CNS-related functions 
(Tullio et al. 2000). Interestingly, CCK-2R2 in brain was 
found to be much higher expressed than the other two 
receptors identified in this study. Further studies required 



to understand the functional distinction (if any) between 
CCK-2R1 and CCK-2R2. 

Perspectives and significance 

In the present study, three CCK-R receptor genes have 
been identified and cloned from brain and gut tissues of 
Atlantic salmon. Within gut segments, CCK-1R mRNA 
were highly expressed in pancreas whereas CCK-2Rl/gas- 
trin expression was more abundant in the midgut. The 
higher levels of expression of CCK-2R2 in brain com- 
pared to other organs, together with the observation that 
it is the most abundant transcript among the three 
CCK-R genes analyzed in this tissue, is of particular 
interest and further studies are required to establish its 
distinct roles in physiological brain functions. Here we 
also described the major structural and functional motifs 
of these receptors and showed that they are conserved 
to that of mammals. Overall, our studies provide a rele- 
vant contribution toward the understanding of receptor- 
mediated execution of CCK-R biological responses in 
Atlantic salmon, and give an essential foundation to 
further elucidate by comparative studies the potential 
conserved roles of these hormone receptors during verte- 
brate evolution. 
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